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Abstract 

Scalar MSSM Higgs boson production via gluon fusion gg ^ h, H is mediated by 
heavy quark and squark loops. The higher order QCD corrections to these processes 
turn out to be large. The full supersymmetric QCD corrections have been calculated 
recently. In the limit of large SUSY masses a conceptual problem appears, i.e. the 
proper treatment of the large gluino mass limit. In this work we will describe the 
consistent decoupling of heavy gluino effects and derive the effective Lagrangian for 
decoupled gluinos. 



1 Introduction 



Higgs boson [1] searches belong to the primary motivations for present and future colhd- 
ers. In the MSSM two isospin Higgs doublets are introduced for the masses of up- and 
down-type fermions [2]. After electroweak symmetry breaking five states are left as ele- 
mentary Higgs particles, two CP-even neutral (scalar) particles h, H, one CP-odd neutral 
(pseudoscalar) particle A and two charged bosons H^. At leading order (LO) the MSSM 
Higgs sector is fixed by two independent input parameters which are usually chosen as the 
pseudoscalar Higgs mass Ma and tg/3 = V2/V1, the ratio of the two vacuum expectation 
values. Being lighter than the Z boson mass at LO, the one-loop and dominant two-loop 
corrections shift the upper bound of the light scalar Higgs mass to M^ < 140 GeV [3] . The 
couplings of the various neutral Higgs bosons to fermions and gauge bosons, normalized 
to the SM Higgs couplings, are listed in Table 1. The angle a denotes the mixing angle of 
the scalar Higgs bosons h, H. An important property of the bottom Yukawa couplings is 
their enhancement for large values of tg/3, while the top Yukawa couplings are suppressed 
for large tg/? [4]. Thus, the top Yukawa couplings play a dominant role at small and 
moderate values of ig(3. 
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Table 1: iiiggs couplings in the MSSM to fermions and gauge bosons [V — W,Z] relative 
to the SM couplings. 



Usually the scalar superpartners fL,R of the left- and right-handed fermion components 
mix with each other. However, in this work we will neglect mixing effects. Thus the masses 
of the sfermion states f^R are simply given by 

Kn = "4 + Kr (1) 

where mj denotes the fermion mass and Mi/ji the soft supersymmetry-breaking sfermion 
mass parameters. The neutral scalar Higgs couplings [H = h, H] to non-mixing sfermions 
read [5] 

9lh = 9lf,-rn}gy 

ah - (2) 
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with the couphngs hsted in Table 1. D terms have been neglected in these expres- 
sions. It is important to note that supersymmetry relates the diagonal couplings to the 
corresponding fermion Yukawa coupling involving the fermion mass m^. 

At hadron colliders as the Tevatron and LHC neutral Higgs bosons are copiously 
produced by the gluon fusion processes gg h/H/A, which are mediated by top and 
bottom quark loops as well as stop and sbottom loops for the scalar Higgs bosons H in 
the MSSM (see Fig. 1) [6]. 




Figure 1: MSSM Higgs boson production via gluon fusion mediated by top- and bottom 
quark as well as stop and sbottom loops at leading order. 

The pure QCD corrections to the (s)top and (s)bottom quark loops are known includ- 
ing the full Higgs and (s) quark mass dependences [7]. They increase the cross sections by 
up to about 100%. The limit of very heavy top quarks and squarks provides an approxima- 
tion within ~ 20 — 30% for tg/? < 5 [8]. In this limit the next-to- leading order (NLO) QCD 
corrections have been calculated [9] and later the next-to-next-to-leading order (NNLO) 
QCD corrections [10]. The NNLO corrections lead to a further moderate increase of the 
cross section by ~ 20 — 30%, so that the dominant part are the NLO contributions. An 
estimate of the next-to-next-to-next-to-leading order (NNNLO) effects has been obtained 
[11] indicating improved perturbative convergence. Moreover, the full SUSY-QCD cor- 
rections have been derived for heavy SUSY particle masses [12-15] and recently including 
the full mass dependence [16]. Ref. [12] addresses the limit of large gluino masses for 
the scalar Higgs couplings to gluons for degenerate squark masses, i.e. without mixing, 
as a special limit of their final result. The result develops a logarithmic singularity for 
large gluino masses which seems to contradict the decoupling of the gluino contributions 
according to the Appelquist-Carazzone theorem [17]. In the pseudoscalar Higgs case this 
logarithmic divergence for large gluino masses does not appear [14]. This work describes 
the resolution of this problem and a consistent derivation of the effective Lagrangian after 
decoupling the gluino contributions. 

The paper is organized as follows. In Section 2 we derive the effective Lagrangian in 
the limit of heavy quark, squark and gluino masses, where the gluinos are much heavier 
than the quarks and squarks in addition. Section 3 summarizes and concludes. 
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2 Decoupling of the Gluinos 



For the derivation of the effective Lagrangian for the scalar Higgs couphngs to gluons we 
will analyze the relation between the quark Yukawa coupling Xq and the Higgs coupling 
to squarks Xq in the limit of large gluino masses in detail. To set up our notation we will 
define these couplings at leading order in the case of vanishing squark mixing as 

2^ . 2,3!^ = nXl 

V V 

(3) 

where v — l/\J\/2Gp ~ 246 GeV denotes the vacuum expectation value of the Higgs 
sector which is related to the Fermi constant Gp. These couplings mediate the scalar 
Higgs decays into heavy quark pairs Ti QQ and into squark pairs Ti — > QQ (see 
Fig. 2a,b). In the following we will derive the modified relation between these couplings 
for scales helowthe gluino mass Mg. This will proceed along the following lines: (i)We will 
start with the unbroken relation between the running MS couplings of Eq. (3) for scales 
above the gluino mass and the corresponding renormalization group equations, (ii) If the 
scales decrease below the gluino mass the gluino decouples from the renormalization group 
equations. This decoupling leads to modified renormalization group equations which differ 
for the two couplings Xq and kXq. This implies that the two couplings deviate for scales 
below the gluino mass, while they are identical for scales above the gluino mass, (iii) The 
proper matching at the gluino mass scale yields a finite threshold contribution for the 
cvohition from the gluino mass to smaller scales, while the logarithmic structure of the 
matching relation is given by the solution of the renormalization group equations below 
the gluino mass scale. We will determine these ingredients in detail in this letter. The 
matching relations can be obtained from the gluino contributions for heavy gluino masses 
in the limit of vanishing external momentum transfers [18] as we will discuss in the next 
sections in detail. 

2.1 ^ QQ 

The gluino contribution to the Higgs vertex with heavy quark pairs is depicted in Fig. 2c. 
We use dimensional regularization in n = 4 — 2e dimensions for the evaluation of the 
one-loop contributions. Since all virtual particles are massive, there are no infrared di- 
vergences. Although dimensional regularization requires the introduction of anomalous 
counter terms in general to restore the supersymmetric relations between corresponding 
couplings and masses, these gluino contributions are free of these terms. The result of the 
vertex contribution in the heavy gluino mass limit vanishes. 
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Figure 2: Scalar MSSM Higgs boson couplings to heavy quarks Q and squarks Q: (a) TiQQ 
coupling at LO, (h) HQQ coupling at LO, (c) gluino contribution to the TiQQ coupling 
at NLO and (d) gluino contribution to the TiQQ coupling at NLO. [H = h,H] 



The full correction to the bare Yukawa coupling requires the addition of the wave function 
renormalization constant Z2 which can be derived from the derivative of the corresponding 
self-energy diagram shown in Fig. 3a. In the limit of large gluino masses the gluino 
contribution is found to be [C^ = 4/3] 

(5) 




The gluino contributions to the bare quark Yukawa coupling can now be derived as 



This result can be used to construct the renormalization group equation for the quark 
Yukawa coupling without the gluino contribution. The full renormalization group equation 
for the running MS coupling including gluon and gluino contributions is at leading order 
given by [19] 



It describes the scale dependence of the running MS coupling for scales larger than the 
quark, squark and gluino masses. However, if the gluino mass is large compared to the 
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Figure 3: Gluino contributions to the self-energies of (a) quarks Q and (h) squarks Q at 
NLO. 



chosen renormalization scale, the gluino has to be decoupled from the renormalization 
group equation. This can be performed consistently by a momentum subtraction of the 
gluino contribution for vanishing momentum transfer while treating the quark and gluon 
contributions in the usual MS scheme [18]. It relates the momentum-subtracted Yukawa 
coupling Xq^mo to the full MS coupling Ag in the following way 




Differentiating this relation with respect to fi'j^ yields the renormalization group evolution 
of the momentum-subtracted coupling corresponding to the low-energy result without the 
gluino, 

dj? ~ ~4 ^ ■^Q,MO[fJ'R) 

Since the matching of the effective theory below the gluino mass scale to the full MSSM 
will be performed at fiR = Mg, Eq. (8) determines the required threshold contribution, 
too, 



Aq,mo(M,) = Aq(M,) |i - ^C'p^^} (10) 

In Ref. [12] the quark Yukawa coupling has been renormalized by introducing the quark 
pole mass rriQ. In the effective theory below the gluino mass scale the running MS coupling 
is related to the quark pole mass as [20] 

9S^ = W'"«){i + c.^j (11) 



2.2 4>^QQ 



The analogous calculation has to be repeated for the Higgs coupling to squarks Xq. The 
corresponding gluino contribution is shown in Fig. 2d. In this case no anomalous counter 



5 



terms are needed, too. The result of the vertex correction in the heavy gluino mass hmit 
is given by 

Here we concentrate only on the diagonal terms, i.e. the 'HqlQl T^QrQr couplings, 
in the no-mixing case since these will be treated by the renormalization of the Yukawa 
coupling. In the general case including squark mixing effects the additional diagonal and 
non-diagonal contributions will be absorbed by the renormalized trilinear coupling Aq. 
The gluino contribution to the squark wave function renormalization constant Z2 for large 
gluino masses can be derived from the derivative of the corresponding self-energy diagram 
as depicted in Fig. 3b, 



a.. 



ATT/i^Vr 1 1 



z,_i^c.^r(i + .)^^j (13) 

In the case of squark mixing the additional non-diagonal contributions arising from the 
self-energy diagram of Fig. 3b will be absorbed by the renormalized mixing angle 9q. As 
for the quark Yukawa coupling the gluino contribution to the Higgs coupling to squarks 
can now be determined. 



For the running MS Higgs coupling to squarks the full renormalization group equation 
including gluon, squark and gluino contributions can be expressed as [19] 

I^R Q^2^ - -'^F Agl^/iRj {lb) 

Note that by virtue of supersymmetry the beta function is twice as large as the corre- 
sponding one for the quark Yukawa coupling in Eq. (7). Eq. (15) describes the evolution 
for scales above the quark, squark and gluino masses. The momentum-subtracted cou- 
pling Aq j^q of the effective theory without gluinos is now related to the full MS coupling 

by 

^Q,MO W = ^{^^R) |l - Cf'^ (]^ log ^ - ^) I (16) 



where the additional contributions of quarks, gluons and squarks to Xq^mo are still MS 
subtracted. The evolution of the momentum-subtracted coupling is fixed by the renor- 
malization group equation, which can be derived by differentiating Eq. (16) with respect 
to 



— di4, — - -^^^ W W (17) 
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This renormalization group equation differs from the corresponding renormahzation group 
equation of the squared momentum-subtracted Yukawa couphng Aq,mo of Eq. (9). Thus 
the supersymmetric relation between these two couphngs is violated for scales below the 
gluino mass where the gluino contribution has to be deleted from the corresponding beta 
functions. This difference is expected, since for heavy decoupled gluinos the residual con- 
tributing particle spectrum is not supersymmetric any more. Moreover, for the matching 
scale = Mg, Eq.(16) determines the threshold correction for the evolution below the 
gluino mass scale, 

Aq,mo(^^) = Aq(M,) |i + ^C^^^} (18) 

Due to the decoupling of the gluinos the soft supersymmetry breaking induces a hard 
supersymmetry breaking at low energy scales [21] as can be inferred from the different 
threshold corrections and the different renormalization group equations below the gluino 
mass scale. 

2.3 Decoupling of gluinos 

Now we are in the position to derive the effective low-energy scalar Higgs coupling to 
gluons. For the consistent decoupling of heavy gluinos their contribution has to be treated 
in the momentum-subtraction scheme as described before. The relation between the 
momentum-subtracted quark Yukawa coupling and the scalar Higgs coupling to squarks 
can be determined as follows. In the supersymmetric theory, i.e. for scales above the 
gluino mass, the supersymmetric relation holds for the running MS couplings, 

i^>^qM^>^qM (19) 

Using Eqs. (8, 16) this leads to a non-supersymmetric relation between the momentum- 
subtracted couplings, 

^>^Q,Moi^Q) = >^Q,MoimQ) [^ + Cf^ (log ^ " ^) } (20) 

In this equation we have set the renormalization scale equal to the quark mass, since at 
this scale the momentum-subtracted Yukawa coupling Xq^mo is related to the quark pole 
mass. Using Eq. (11) for the relation between the quark pole mass and the MS Yukawa 
coupling of the effective theory below the gluino mass scale we arrive at the relation 

= ^«,Mo(™«) {^ + CA (log ^ + i) ) (21) 

The proper scale choice for the effective Higgs coupling to squarks, however, is the squark 
mass. This choice is relevant for an additional large gap between the quark and squark 
masses. Using the renormalization group equation of Eq. (17) we end up with the final 
relation 
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The radiative corrections to the relation between the effective couplings after decoupling 
the gluinos modify the final result of Ref. [12]. This modification can be discussed in 
terms of the effective Lagrangian in the limit of heavy squarks and quarks, 
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where = vXq f^Q{mQ) /rn^. In Ref. [12] the leading term of the supersymmetric 
coefficient for large gluino masses has been derived for equal squark masses ttiq = Mq^ = 
Mq^ as 

C^' -^^-So^^^ 210^""^ (241 

where the mismatch between the couplings of Eq. (22) for scales below the gluino mass 
has not been taken into account, i.e. keeping the relation — 2g^mQ/rn?Q after renor- 

malization as in the supersymmetric limit. Expressing g^ in terms of ^Q^Moi^q) instead, 
this mismatch leads to the additional contribution 

AC5Qcz. = T.log-f + 21og-f + - (25) 
o rriQ rriq o 

Adding both contributions Csqcd = C^qq^j + ACsqcd determines the supersymmetric 
contribution to the effective Lagrangian, 

Csqcd = -g- (26) 

The resulting effective Lagrangian is well-defined in the limit of large gluino masses and 
thus fulfills the constraints of the Appelquist-Carazzone decoupling theorem [17]. 

The resulting effective Higgs coupling mo^'^^q) ^^"^ determined by solving the 
renormalization group equations Eqs. (9, 17) which are valid for scales below the gluino 
mass. Taking into account the proper matching to the fully supersymmetric MS coupling 
at the scale jiR — Mg according to Eqs. (10, 18) we arrive at the expression 




'^Q,MO\"''Q 



h.Mo(^6) = 2gg^ ^fr„,,- , (27) 



where Pq — (33 — 2Nf — Np)/V2 denotes the leading order beta function of the strong 
coupling Us- Np is the number of contributing quarks and Np the number of con- 
tributing squark flavours. Decoupling only the gluino from the supersymmetric spectrum 
[Np = Np = 6] its value is given by /Sq = 5/4. The expression (27) can be used to evaluate 
the effective coupling Xq Moi^g) i^^om the quark pole mass to leading logarithmic accu- 
racy. This expression resums the leading logarithms of the gluino mass and provides the 
consistent matching of the low-energy Lagrangian to the full MSSM. In this way the value 
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of the gluino mass is still measurable as the scale at which the two couplings Xq j^q and 
kAq merge after taking into account the corresponding threshold corrections. Note 
that the pure perturbative expansion of Eq. (27) reproduces Eq. (22) up to 0{as). Finally 
we would like to emphasize that Eqs. (22, 27) do also hold in the general mixing case for 
the part of the Higgs couplings to squarks which is related to the quark Yukawa coupling, 
if the individual squark masses are chosen as the renormalization scales. 

Since pseudoscalar Higgs bosons only couple to different squark states, i.e. to QlQr 
and vice versa, there are no squark loops at LO so that the pseudoscalar coupling to 
squarks contributes at NLO for the first time. This explains, why the logarithmic singu- 
larity for large gluino masses does not appear in the pseudoscalar case at NLO [14]. 

3 Conclusions 

In this work we have described the consistent decoupling of heavy gluino mass effects from 
the effective Lagrangian for the scalar Higgs couplings to gluons within the MSSM. We 
have shown that a careful extraction of the gluino contributions in the large gluino mass 
limit induces a modification of the supersymmetric relations between the Higgs couplings 
to quarks and squarks. This modification involves non-trivial logarithmic gluino mass 
contributions to the effective couphngs. They exactly cancel the left over gluino mass 
logarithms of the previous work of Ref. [12] which did not take into account the mismatch 
between the Higgs couplings to quarks and squarks at scales below the gluino mass. This 
work ensures that the gluino contributions decouple explicitly for large gluino masses in 
accordance with the Appelquist-Carazzone theorem [17]. The gluino mass remains in the 
effective low-energy theory as the matching scale to the full MSSM. Analogous methods 
have to be applied to all processes if decoupling properties are analyzed. 
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